Nonequilibrium Landau-Zener-Stiickelberg spectroscopy in a double quantum dot 
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We study theoretically nonequilibrium Landau-Zener-Stiickelberg (LZS) dynamics in a driven 
double quantum dot (DQD) including dephasing and, importantly, energy relaxation due to en- 
vironmental fluctuations. We derive effective nonequilibrium Bloch equations. These allow us to 
identify clear signatures for LZS oscilations observed but not recognized as such in experiments [Pe- 
tersson et al., Phys. Rev. Lett. 105, 246804, 2010] and to identify the full environmental fluctuation 
spectra acting on a DQD given experimental data as in [Petersson et al., Phys. Rev. Lett. 105, 
246804, 2010]. Herein we find that super-Ohmic fluctuations, typically due to phonons, are the main 
relaxation channel for a detuned DQD whereas Ohmic fluctuations dominate at zero detuning. 
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Quantum electronic devices, as qubits realized by dou- 
ble quantum dots (DQD), require coherence times which 
exceed their quantum operation time during which the 
DQD is typically strongly driven by external voltage 
pulses. Tremendous research efforts studied semiconduc- 
tor based devices to achieve coherent quantum control [l|- 
0] • Many fluctuation sources of the noisy solid state envi- 
ronment, which act on the electron in the DQD and thus 
destroy coherence, were revealed but a comphrchensive 
picture is elusive. Furthermore, driving by voltage pulses 
causes an intrinsic nonequilibrium situation in which re- 
laxation competes with driving inn which renders a 
theoretical description of the dissipative nonequilibrium 
dynamics highly nontrivial. 

Here, we theoretically study the dissipative nonequi- 
librium dynamics of a single electron charge qubit de- 
fined in a DQD embedded in a noisy solid state envi- 
ronment driven by voltage pulses. While DQD charge 
qubits have relatively short coherence times, this disad- 
vantage is compensated by the possibility of fast quan- 
tum operations. We model the DQD and its dissipa- 
tion as a quantum two-level system in an open quan- 
tum system approach (llT |. We determine the dissipa- 
tive nonequilibrium real time dynamics (initialized by 
applying voltage pulses) by deriving effective nonequilib- 
rium Bloch equations (NBEs). These allow fast numer- 
ical treatment in contrast to numerical exact methods 
Q and thus allow a comphrehensive analysis of re- 
cent experiments by Petersson et al. [l[ and Dovzhenko 



et al. 12l | . In these experiments quantum control of a 



single electron was achieved by means of applying ultra 
short voltage pulses to control gates of the laterally de- 
fined DQD. In these time ensemble measurements the 
DQD was cycled (with 40 MHz repetition rate) between 
two different ground state configurations while its aver- 
age charge occupation was continuously detected via the 
electric current through a capacitively coupled quantum 
point contact (QPC). The applied voltage pulses gener- 



ate Landau-Zener-Stiickelberg (LZS) dynamics [l3| and 
we can identify so far unexplained features in the exper- 
imental data as signatures of coherent LZS oscillations. 

In the experimental ensemble measurements the de- 
phasing time due to slow noise 0, [H| is much shorter 
than relaxation times T\ and thus dominates the decoher- 
ence times T 2 since T 2 _1 = + (2T 1 )~ 1 . Relaxation 

and dephasing could be caused by thermal phonons, in- 
trinsic or externally triggered charge noise flBl Hit , or 
detector back-action 



12| neglected re 



[16H22J]. Ref. [1, _ 
laxation and employed solely an 1// noise model 23l- 27j 
for dephasing which allowed them to describe the detun- 
ing dependence of the observed decoherence of the DQD. 

The experimentally observed steady state charge occu- 
pation of the DQD is, however, heavily influenced by re- 
laxation [l(3|. Therefore, in this letter, we go beyond this 
simple dephasing model and include various environmen- 
tal fluctuation spectra to describe dephasing and relax- 
ation. This allows us to simulate the full nonequilibrium 
real time dynamics experimentally studied. We find that 
the observed visibilty reduction is caused by relaxation. 
Moreover, by analyzing the measured LZS dynamics we 
identify the full environmental fluctuation spectrum as a 
sum of three processes. In addition to slow noise, already 
considered in ref. which causes detuning dependent 
dephasing, super-Ohmic fluctuations, as typically origi- 
nating from phonons, are the main relaxation channel for 
the detuned DQD. Near zero detuning, however, Ohmic 
fluctuations dominate relaxation. The latter also limit 
the experimentally observed maximal decoherence time 
of T 2 - 7 ns. 

Modelling of pulse driven DQD We model the single 
electron DQD using a two-level system Hamiltonian 



H 



h{t)a z 



(1) 



where the eigenstates of <j z correspond to the electron 
in the left / right dot, e(t) is the level detuning and 
A is the interdot tunnel splitting [28|. Experimentally, 




FIG. 1. (Left): Energy detuning profile for a voltage pulse 
with rise time t r , pulse plateau time t p , pulse duration t v — 
t p +2t r - If cycled the pulse periodically repeats with repetition 
time t Icp . (Right): Decoherence rate versus detuning e p . Full 
theory (black full line), Ohmic dephasing (red dashed line) 
and decoherence due to relaxation (blue dot-dashed line). 

gate voltage pulses are applied to change the level detun- 
ing as sketched in Fig. [TJlcft). Initially the detuning is 
e = e(0) ^> A and the DQD in the according ground 
state (0,1) [electron is in the right dot]. A voltage pulse 
then drives the system within a rise time t r to a plateau 
detuning e p close to resonance (e = 0) and keeps it there 
for a plateau time t p . Then the DQD is driven back 
within t r to the initial detuning and the probability P(i,o) 
of occupation of the excited state (1,0) [electron in the 
left dot] is studied as a function of the pulse duration 
t v = 2t r + t p . For a quantum mechanical two- level sys- 
tem without dissipation we expect coherent oscillations 
in the voltage pulse time t v . 

A voltage change on a single gate P, [l2| not only af- 
fects the level detuning but also causes a common energy 
shift of both states (1,0) and (0,1) which is not relevant 
in the following and thus neglected. The sharp corners 
of e(t) in Fig. QJleft) are smoother in reality and might 
even contain oscillatory features (due to a finite band- 
width transfer function). At voltage pulse times t v > 2t r 
deviations between linear voltage ramps and more accu- 
rate descriptions are negligible. 

Solving the quantum dynamics of the DQD with a sin- 
gle applied voltage pulse (details given below) results in 
the probability Pnfi) (t — t v ) after the voltage pulse, plot- 
ted in Fig. [2] as a function of t v and plateau detuning 
t p . In order to model the experiments 0, [l2| we use 
A/h = 4.5 GHz, e = 200,ucV ~ 11A and t r = 35 ps 
which cor resp onds to the fastest experimental achievable 
rise time [29[. The sweep speed v p = \e p — eo\/t r can be 
approximated by Vq = |eo|/i r for eo e p [slope in Fig. 
QJleft)] which provides an indication on the overall adia- 
baticity of the dynamics [Hj]. For v = 200//eV/35ps ~ 
llA 2 /fi,, Fig.[2]shows coherent oscillations between states 
(0,1) and (1,0) with frequency E/h and eigenencrgy 

E = ^A 2 + e 2 of Eq. ([J). At zero detuning P(i, )(t v ) 

oscillates between and 1. With increased e„ the os- 



2 




500 1000 1500 2000 

t v = t p +2t r [ps] 

FIG. 2. P im (t„) for the fully coherent case (Ti = T 2 = T 2 * = 
oo), right after a single voltage pulse, vs. level detuning e p and 
pulse duration t v = t p + 2t r for eo = 200/ieV (eo = 2000^ieV) 
for the main figure (inset). 

cillation frequency grows while its amplitude gradually 
decreases. Fig. [5] shows a clear asymmetry: the visibility 
is larger for e p < compared to e p > 0. Remarkably, this 
asymmetry was observed experimentally, too (but not 
explained) [l| . The asymmetry decreases with increasing 
sweep speed as highlighted by the inset of Fig. [2] which 
shows an almost symmetric -P(i,o)(^) f° r e o = 2000 /xeV 
leading to v > 110A 2 /?i. We conclude that the asym- 
metry is solely an effect of adiabaticity when the DQD 
is driven through the avoided crossing at zero detuning. 
For e p < 0, the quantum system accumulates, during the 
voltage pulse, not only phase due to the coherent oscilla- 
tion of the electron between the two dots but also due to 
the superposition state occupied between two successive 
Landau-Zener transitions [13| . As such, the asymmetry 
is a clear signature of coherent LZS oscillations. 

Periodically Cycled Pulses In an ensemble measure- 
ment with continuous charge detection, as in Ref. [lj with 
t rep = 25 ns, it is essential to choose i rep ;§> t v to ensure 
readout of the charge occupation after application of the 
pulses (of duration t v < T2). An interpretation in terms 
of an ensemble measurement, which simply averages over 
many independent shots, further requires t rop 3> I2 and, 
interestingly, i rep ~ Ti, where Ti is the (thermal) en- 
ergy relaxation time which depends on the detuning: for 
^rcp ^ Ti initialization into configuration (1,0) is guar- 
anteed, but the visibility of the continuous measurement 
is close to zero as mostly (1,0) is occupied; for t rop <C Ti 
initialization into (1,0) independently of -P(i,o) -right af- 
ter the pulse- is impossible. No matter of the choice of 
pulse sequence, continuously cycled pulses will result in a 
steady state which, in principle, contains the information 
of dephasing and energy relaxation times [l(J • 

In order to model the experimental repeated pulse 
train sequence we repeat our simulation after the first 
cycle (up to the repetition time t rcp ) with the final statis- 
tical operator of the previous cycle as initial state. This 
procedure is repeated until the population -P(i,o) [mca- 
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t v = t p +2t r [ps] 

FIG. 3. (Main): P(i.o) [(Inset): P(i,o)(tv) with color scale 
from to 1] versus detuning e p and voltage pulse time 
t v = t p + 2t r for a driven DQD with eo = 200^ieV includ- 
ing relaxation and dephasing. 

sured as average over a full cycle as plotted in Fig.fljlcft)] 
changes by less than 0.005. This -P(i,o) approximates the 
experimentally observed steady state population. 

Driven dissipative dynamics Including dissipation 
into the driven dynamics of the DQD causes a compe- 
tition between driving and relaxation [1, [j| which ren- 
ders all but expansive numerical treatments inadequate. 
Following the standard approach within open quantum 
dynamics (llj we couple the driven two-level Hamiltonian 
([T]) of the DQD to environmental fluctuations described 
by harmonic oscillators. This results in 

ff to t =H(t) + ^-J2 X k(Pk + bl) + Y, UJ kb\b k (2) 

k k 

with bosonic annihilation/creation operators bk/bj^. The 
spectrum Gtui) = ^ fc Xf,S(uj — ujf,) is typically a smooth 
function at the energies of interest, i. e. G(uj) ~ 
2aui s exp(— l/j/lo c ) with spectral exponent s, cut-off fre- 
quency uj c and the coupling strength a. Typically [l]], 
A ^ 72 S> 71 with the relaxation rate 71 = 1/T\ 
and the dccohercncc rate 72 = I/T2, which justifies 
a weak DQD-cnvironmcnt coupling approach. Due to 
the time dependence of H{t) the weak coupling Born- 
Markov approximation, however, fails. Assuming that 
the change of energy during the memory time in the en- 
vironment is small, allows an additional adiabatic rate 
approximation. For this, we switch to the time depen- 
dent basis where H (t) = ^E(t)r x is diagonal with Pauli 
matrizes Tj describing the eigenstates with energy dif- 
ference E(t) = y/ A 2 + e(t) 2 . Employing a lowest or- 
der Born approximation for the memory kernel in the 
system- bath coupling [3(1 (3l| , one obtains for the time 
evolution of the components of the statistical operator, 
pit) = i(ll + r j T j)i noncquilibrium Bloch equations 

d t r x (t) = <P'{t)r z {t) -TlW^W-rf (*)] 
d t r y (t) = -j 2 (t)r v (t) -E(t)r z (t) (3) 
d t r z {t) = E(t)r y (t) -j 2 (t)r z (t) -<j>'(t)r x (t) 



A/h = 3.3 GHz - A/h = 4.5GHz -- A/h = 6.6 GHz 
f s /\ /. 




FIG. 4. P(i,o) for zero detuning versus voltage pulse time 
t v = t p + 2t r for various tunnel couplings. 

with time dependent momentary equilibrium r x q (t) = 
tanh[/3 \E{t)} and time dependent decay rates 

7i W = 7fc coth[\f3E(t)]u 2 {t)G[E{t)\ (4) 

I2(t) = \n{t) +T 2 (t) (5) 

where T 2 {t) = 1/T£(t) is the dephasing time. Herein, 
4>{t) = arctan[e(i)/A], u(t) = cos 4>(t), and v(t) = 
sin (j)(t) . The presented nonequilibrium Bloch equations 
with time-dependent rates and momentary equilibrium 
describe adequately dissipative Landau- Zener dynamics, 
i.e. the regime of competition between driving and relax- 
ation, for weak DQD-cnvironment coupling. We ensured 
this by extensive comparisson with numerical exact re- 
sults [H, IH . Details will be presented elsewhere. 

To proceed, we need the spectra of fluctuations acting 
on the DQD. Estimates can be gained by the observations 
for fixed DQD parameters of Petersson et al. Q of a 
relaxation time Ti(eo) = 10 ns at the initial detuning eo, 
a decoherence time of T 2 ~ 7 ns at zero detuning and the 
detuning dependent dephasing times presented in Fig. 4c 
of Ref. [l|; all at the temperature T = 80 mK. 

Relaxation rates Charge fluctuations and bulk 
phonons likewise couple to the DQD. In order to include 
both, we consider an Ohmic fluctuation spectrum with 
s = 1 (charge fluctuations) as well as a super-Ohmic one 
with s = 3 (phonons) The rate (@| reflects a one- 

boson process and accordingly only yields a substantial 
relaxation if the spectrum at the eigenenergy of the DQD 
is finite, i. e. E(t) <C ui c . 

At fixed T, A and eo we can determine the coupling 
strength a (assuming E(t) -C lo c ) using the known energy 
relaxation time Ti(eo), and estimate Ti(e = 0) using Eq. 
dU) and u(t ) = A/E (t): For Ohmic fluctuations Tf =1 (e = 
0) = [A/^/A 2 + el) ■ Ti(e ) while for super-Ohmic fluc- 
tuations rf =3 (e = 0) = (VA 2 + eg/A)-Ti(e ). With the 
reported Ti(e ) = 10 ns we find T[ =1 {e = 0) ~ 900 ps as- 
suming Ohmic fluctuations compared to T-f =3 (e = 0) ~ 
110 ns assuming super-Ohmic fluctuations. According 
to Eq. (|5|) energy relaxation causes an upper bound of 
the decoherence T 2 < 2T\. The decoherence time of 
T 2 = 7 ns, observed in the experiment at e ~ [H, lies 
well in between our predictions T( =1 <C 7ns <C T-f =3 . 
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An Ohmic fluctuation spectrum alone -consistent with 
Ti(eo) = 10 ns- would result in much too fast decoher- 
ence at e = compared to the experimental results. 
Hence, super-Ohmic fluctuations, typically caused by 
phonons, are the main relaxation mechanism at eo- 

Decoherence and Dephasing rates For a super-Ohmic 
spectrum T 2 (t) = 0, and the decoherence rate 72 [see 
Eq. ([5])] would be solely determined by the energy relax- 
ation mechanism. The spectrum of charge noise, which 
is well known to cause additional dephasing, strongly de- 
pends on the sample and how it has been treated [lij ]. 
Often, charge noise can be assumed to be slow noise as 
done by Pctcrsson et al. [lj. They use an 1/ f dephas- 
ing model [25f which typically results from background 
charge noise [23, 12411 and might be described in terms of 
sub-Ohmic noise |2(j |27| with s = and uj c <C A. Such 
slow noise is a major dephasing source in realistic devices 
and causes solely dephasing. We aim at unraveling the 
relaxation mechanisms present rather than the dephasing 
sources. Relaxation is not influenced by slow noise and 
thus we describe slow noise here using a simplified Ohmic 
slow noise model resulting in T 2 (t) = (4ir / 'h)v 2 '(t)a 2 k-QT 
characterized by a coupling strength 02 which we fit 
to the experiment. Our Ohmic slow noise model cap- 
tures two important points, namely, that it does not 
influence relaxation, i.e. u> c2 <C A, and that it cou- 
ples to the position of the charge in the DQD leading 
to r 2 (t) cx v 2 {t) = e 2 (i)/(A 2 + e 2 (i)), which coincides 
with the dependence used by Petersson et al. jl| in their 
l/f noise model. As typically done, we, herein, neglect 
fluctuations in A. 

The red dashed line in Fig.QJright) represents the con- 
tribution to decoherence by our slow Ohmic dephasing 
noise and thus reflects most of the decoherence (black 
full line in Fig. QJright) which reproduce the measured 
data in Fig. 4c of Ref. [l( remarkebly well) except at very 
small dctunings e p A. For v(t) = e p = 0, however, the 
slow noise dephasing rate T 2 (t) = and 72 (t) = ^71 (t). 
Super-Ohmic fluctuations result in T-f =3 (e = 0) ~ 110 ns, 
much longer than the zero detuning decoherence time 
T-2 = 7 ns actually measured. To resolve this discrep- 
ancy, we add a third noise source, namely weak Ohmic 
fluctuations which contribute to relaxation at zero de- 
tuning in addition to the super-Ohmic contribution of 
phonons. The contribution of weak Ohmic fluctuations 
to energy relaxation at large detuning is very small. The 
physical origin of this mechanism could be related to fast 
local potential fluctuations (e. g. via voltage noise on the 
lead gates). 

Our complete model includes three fluctuations spectra 

3 

G(u) = 2a 1 uj 3 e- ul/u] -- 1 + ^ lu^e-"'^ (6) 

with u! c i,uj C 3 ^ A and uj c2 <C A. Choosing aiA 2 = 
8.09 • 10~ 5 , a 2 = 3.53 • 10" 2 and a 3 = 2.73 ■ 10~ 3 for 



the coupling strengths reproduces the experimental ob- 
servations [l( of Ti(eo) = 10 ns, T 2 (e = 0) ~ 7ns and the 
decoherence time versus detuning with measured data 
in Fig. 4c of Ref. [l| and our predicted decoherence rate 
plotted as black full line in Fig. QJright). 

Dynamics of the driven dissipative DQD Solving the 
nonequilibrium Bloch equations ([3]) with the fluctuation 
spectrum ((S]) numerically results in Fig. (3[ the main re- 
sult, where we plot the steady state occupation P(i : o) ver- 
sus voltage pulse time and plateau detuning for a tunnel 
coupling A/h = 4.5GHz. In comparison, the inset shows 
the instantaneous probability P(i,o) = with the color 
code stretched by a factor of two, i.e. yellow = 1. The 
continuous current measurement in the QPC reduces the 
visibility for the coherent dynamics by a factor of 2. In 
comparison to the undamped case in Fig. [2] we observe 
an overall reduction in oscillation amplitude. On top, 
the LZS oscillations smear out with both, more negative 
detunings e p and longer voltage pulse times t v : P(i,o) is 
flat in the upper right corner of Fig. [3J The same be- 
haviour was experimentally observed, but not explained, 
by Petersson et al. Our results indicate that the 

LZS oscillations smear out as a result of relaxation but 
the remaining asymmetry with respect to e p is a result 
of adiabaticity. Thus, the experimentally observed fea- 
ture is a clear signature of LZS oscillations. The overall 
qualitative agreement between our simulation results and 
the experimental data [l| is very good. The quantitative 
overestimation of the oscillation amplitude being a factor 
of 2 is likely due to uncertainties regarding experimental 
details [32|. 

Interdot tunnel coupling Petersson et al. [l[ ob- 
serve an increased (decreased) oscillation amplitude for 
a smaller (larger) interdot tunnel coupling, i.e. A/h = 
3.3 GHz (6.6 GHz). We find the same tendency. It results 
from the fact that the relaxation rate [eq. ([4])] increases 
strongly with the tunnel coupling. Fig. [4] plots our pre- 
dictions for P(i,o) as a function of t v for zero detuning for 
the three different tunnel couplings. 

Conclusions We studied theoretically the dissipative 
nonequilibrium dynamics of a single electron DQD driven 
by voltage pulses. We couple the DQD additionally to 
environmental fluctuations causing relaxation and de- 
phasing. Extending standard Born-Markov approaches 
to driven systems, we derive nonequilibrium Bloch equa- 
tions exhibiting time dependent rates and the momen- 
tary equilibrium. This approach allows efficient numer- 
ical simulations of the full nonequilibrium real time dy- 
namics and a comprehensive analysis of the experimental 
data We identify an asymmetric occupation of the 
left dot in respect to the detuning between the dots in 
ref. [l[ as a clear experimental signature of LZS dynam- 
ics. A full analysis of the LZS dynamics furthermore 
allows us to specify the full environmental fluctuation 
spectrum acting on the DQD studied by Petersson et al. 
[l| as sum of three processes. Besides slow noise causing 
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the detuning dependent strong dcphasing, supcr-Ohmic 
fluctuations, as typically originating from phonons, arc 
the main relaxation channel for a detuned DQD. At zero 
detuning, however, Ohmic fluctuations, which might be 
caused by gate voltage noise, dominate relaxation and are 
also the main cause for the decoherence at zero detuning. 
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